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1 Introduction 

The determination of quark masses has been transformed in the past few years by ac- 
curate results from realistic lattice QCD. This has meant a range of new methods for 
both b and c quarks [1] which I will describe along with results from continuum tech- 
niques. The recent improvement in Lattice QCD actions along with the generation of 
gluon configurations including sea quarks on fine lattices has given us viable methods 
for handling heavy quarks with relativistic actions. This has allowed us for the first 
time to connect the heavy and light sectors through accurate determination of ratios 
of quark masses, such as m c J m s [2 J . Such ratios enable us to leverage the accuracy of 
c and b mass determinations into accurate s and light quark masses. The comparison 
of quark masses from different methods and formalisms provides a strong test of QCD 
and the masses themselves are needed in calculations for the cross-sections of various 
processes at LHC, for example H — > bb. 

2 Lattice QCD methods 

Lattice QCD calculations have the advantage over the real world of having direct 
access to the parameters of the QCD Lagrangian. These are tuned against exper- 
imentally well-determined hadron masses, one for each quark mass. The masses in 
the lattice QCD Lagrangian are in a lattice regularisation scheme (which depends 
on the discretisation of the QCD Lagrangian used) and at a scale which depends on 
the inverse lattice spacing. To convert them to the standard MS scheme at a fixed 
scale must be done either directly, by determining a renormalisation constant for the 
conversion, or indirectly, by calculating on the lattice some ultraviolet-finite quantity 
whose value is known in the continuum in terms of the MS quark mass. 



The direct conversion, m(/i) = Z (fia)mi att (l / a) , either calculates Z in lattice QCD 
perturbation theory or matches lattice calculations to a MOM-scheme definition that 
can then be converted to MS using continuum QCD perturbation theory. 

A well-developed and accurate indirect method for heavy quarks is the current- 
current correlator method [3] . Here the continuum limit of time-moments of a well- 
defined lattice QCD heavyonium correlator is compared to continuum QCD pertur- 
bation theory for Q 2 -derivative moments of the heavy quark vacuum polarization 
function. The perturbation theory is known through for some moments since this 
method has been developed for accurate calculations in the continuum - see, for ex- 
ample Pj. There m& and m c are extracted from experimental data by isolating the b 
and c quark contributions to cr(e + e~ — > hadrons) 

Handling b and c quarks in lattice QCD is made difficult by the presence of po- 
tentially large discretisation errors when the heavy quark mass in lattice units, mqa, 
is 0(1). These errors can be beaten down by systematically improving the discretisa- 
tion of the Dirac Lagrangian and this is the best approach for c quarks. For example, 
the Highly Improved Staggered Quark (HISQ) action [5J gives only few % errors for 
c physics at lattice spacings, a « 0.1 fm. Using this same approach for b physics [6] 
requires very fine lattices and extrapolation to the b quark mass from smaller masses 
where discretisation errors are under better control, but has the advantage that the 
same lattice action can be used for all 5 quarks from u to b. Alternatively for b 
quarks there are a number of discretised effective theory approaches that avoid large 
discretisation errors because the b quarks are nonrelativistic. The price to be paid, 
however, is a more complicated action that needs more renormalisation. 
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Figure 1: Comparison plots for (from left to right): m c , m c /m s and m&. Masses are 
in the MS scheme with a scale equal to the mass itself. The ratio is scale- invariant. 
The lattice QCD results marked as 'w, d sea' do not include s quarks in the sea and 
therefore may not be comparable to the other numbers. 
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3 Results 



Figure [T] compares results from lattice QCD and continuum determinations for m c 
(left) and m& (right). For m c 1% errors are possible and for m&, 0.5%. The most 
accurate values are: m c (m c ) = 1.273(6) GeV j6] and ra&(mi,) = 4.163(16) GeV [I]. 
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Figure 2: The figure on the left shows the ratio m c /m s plotted against the square 
of the lattice spacing for results using HISQ quarks on rif = 2 + 1 configurations [2] 
(blue open squares) and Wilson-type quarks on rif = 2 configurations [15] (red open 
triangles). Errors are only shown on the final extrapolated result (filled symbol) at a = 
0. Green bursts show preliminary HISQ results onn/ = 2 + l + l configurations [TB"] . 
The figure on the right shows the ratio m h /m c , divided by the ratio of the pseudoscalar 
meson masses, plotted as a function of the pseudoscalar heavyonium mass in the range 
from nih = m c to rrih = [6]. 

In the plot for m c , the result from the ETM Collaboration |7] is a direct determi- 
nation using the twisted mass quark formalism [5] on gluon configurations including 
only u/d quarks in the sea {rif = 2) and a Z factor determined using the MOM 
scheme. Note that this result is m c {rif = 2) and so may differ from the other results 
since there is no perturbative way of correcting for missing sea s quarks. The other 
results use the current-current correlator method. The HPQCD result [H] uses HISQ 
quarks on gluon configurations including sea u, d and s quarks {rif = 2 + 1) com- 
paring time-moments of the absolutely normalised (in this formalism) pseudoscalar 
charmonium correlators to continuum QCD perturbation theory. The rif = 3 quark 
mass is converted to rif = 4 perturbatively. The lower two results [U |9] are from com- 
paring continuum QCD perturbation theory to moments of the c quark contribution 
(effectively via the vector charmonium correlator) to a(e + e~ — > hadrons). The two 
results differ in their treatment of errors, particularly the perturbative error which is 
handled in both cases by varying the scale at which a s is determined. Both have a 
significant error from the value of a s . The lattice QCD result from HPQCD [6] has 
an advantage here in that a s is simultaneously determined along with m c (and m^). 
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A particularly accurate extraction of a s is possible from the lowest moment of the 
pseudoscalar correlator, which has no explicit mass dependence. In fact the lattice 
calculation fits the 4 lowest moments simultaneously and for a range of masses from 
c to b. This allows both a constraint on higher order terms in perturbation theory 
beyond and a determination of the perturbative error from the effect of including 
them. Moments of the lattice QCD vector correlator can be compared directly to 
the results extracted from cr(e + e~ — > hadrons) [I] and this provides a 1.5% test of 
QCD [TO] (see also [II]). 

The rightmost plot of Figure [T] compares results for m^. The top value from 
HPQCD/HISQ [6] uses the current-current correlator method and is determined in 
the same analysis as m c above. The bottom result is also from the same continuum 
determination as the m c above |4J, using e + e~ data in a different region of \fs. These 
are the most accurate results. The two central lattice QCD results use different 
techniques and so far both include only u/d quarks in the sea, so the results may 
not be comparable. The Alpha collaboration p2] uses static (infinite mass) quarks 
but with 1/rriQ corrections and determines the b mass fom the binding energy of a 
heavy-light meson using nonperturbative step-scaling to determine the energy offset. 
The ETM collaboration [13] interpolates the ratio of the heavy-light meson mass to 
the quark mass between relativistic twisted mass quarks and static quarks. A direct 
determination of mj is available from lattice NRQCD on gluon configurations that 
include u, d, s and c quarks in the sea, but it uses a Z factor only determined to 
0(a s ) in lattice QCD perturbation theory so is rather inaccurate [T4"] . 

Figure [2] shows the quark mass ratios m c /m s and mi,/m c determined from lattice 
QCD calculations. The ratio of lattice QCD masses translates directly to the ratio of 
MS masses (at the same scale) because the Z factors cancel when the same formalism 
is used for both quarks. This is then a fully nonperturbative determination and can 
be done with an accuracy of 1%. 

The leftmost plot shows results for m c /m s as a function of a 2 for HISQ |2j and 
Wilson-type quarks [15J. The improved a-dependence of HISQ quarks is evident, 
particularly for new results onn/ = 2 + 1 + 1 configurations that include a further 
improved gluon action and HISQ quarks in the sea (TBJ. A comparison of different 
lattice QCD results for m c /m s is given in the middle plot of Figure [TJ compared to 
the ratio obtained from the Particle Data Tables [17] running m c and m s to the same 
scale. The best, and only n/ = 2 + 1, result is m c /m s = 11.85(16) from HPQCD [2]. 
It can be used with an accurate value for m c to obtain a 1% accurate m s [5J. 

The rightmost plot of Figure [2] shows HISQ results for md/m c over the range of 
heavy quark masses from c to b. Very little a-dependence is seen and a 1% accurate 
result at the b can be obtained: rrib/m c = 4.49(4) [6j. This agrees well with the ratio 
of masses obtained from the current-current correlator method and provides a good 
nonperturbative test of that method. 
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